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Abstract 


We extend investigations of infrared dynamics in accelerating universes. In the 
presence of massless and minimally coupled scalar fields, physical quantities may ac¬ 
quire growing time dependences through quantum fluctuations at super-horizon scales. 
From a semiclassical viewpoint, it was proposed that such infrared effects are described 
by a Langevin equation. In de Sitter space, the stochastic approach has been proved to 
be equivalent to resummation of the growing time dependences at the leading power. 
In this paper, we make the resummation derivation of the Langevin equation in a gen¬ 
eral accelerating universe. We hrst consider an accelerating universe whose slow-roll 
parameter is constant, and then extend the background as the slow-roll parameter be¬ 
comes time dependent. The resulting Langevin equation contains a white noise term 
and the coefficient of each term is modified by the slow-roll parameter. Furthermore we 
find that the semiclassical description of the scalar fields leads to the same stochastic 
equation as far as we adopt an appropriate time coordinate. 
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1 Introduction 


In accelerating universes, the degrees of freedom at super-horizon scales increase with time. 
This increase gives rise to a growing time dependence to the propagator of a massless and 
minimally coupled scalar held. In some held theoretic models on the backgrounds, physical 
quantities acquire secular growths through the propagator. We call them quantum infrared 
(IR) ehects in accelerating universes. 

For instance, in de Sitter (dS) space, the propagator of a massless and minimally coupled 
scalar held depends logarithmically on the scale factor of the Universe M. By employing 
the Schwinger-Keldysh formalism SIS], we can perturbatively investigate interacting held 
theories on a time dependent background. At each loop level, the IR ehects in dS space 
manifest as polynomials in the logarithm of the scale factor. At late times, the leading 
IR ehects come from the leading logarithms at each loop level. This fact indicates that 
perturbation theories eventually break down if the dS expansion continues long enough. 

From a semiclassical description of the scalar held, A. A. Starobinsky and J. Yokoyama 
proposed that the IR ehects can be evaluated nonperturbatively by a Langevin equation m- 
N. C. Tsamis and R. P. Woodard proved that the stochastic approach is equivalent to 
the resummation of the leading logarithms [8]. These investigations were performed in 
interacting held theories with potentials in dS space. 

In a general accelerating universe, it holds true that the propagator of a massless and min¬ 
imally coupled scalar held is growing with time. Therefore, in a similar way to dS space, 
some nonperturbative analyses of the IR ehects are necessary in interacting held theories 
with the scalar held. 

In this paper, we extend the resummation formula of the leading IR ehects in a general 
accelerating universe. Although there have been some attempts to derive the stochastic 
equation in a general accelerating universe ISM, the consistency with the resummation 
formula has not been verihed completely. We conhrm that in a general accelerating universe, 
the Yang-Feldman formalism is reduced to a Langevin equation with a white noise term up 
to the leading IR ehects. Furthermore the resulting Langevin equation shows that in the 
semiclassical description of the scalar held, the choice of time coordinate is crucial to obtain 
the correct coefficient of its nonlinear terml§- 

As a hrst step, we consider an accelerating universe whose slow-roll parameter is constant. 
On the background, the wave function of a massless scalar held is expressed by the Bessel 
function. As utilized in [T3HT6], this fact facilitates analytic investigations. In Section [21 we 
review the free scalar held theory on the background, especially its IR dynamics. In Section 
El we hrst derive the Langevin equation in the accelerating universe by the resummation of 

^ The dependence of the choice of time coordinate is discussed for the stochastic dynamics of inflatons m- 
Our claim is that it is crucial even for the stochastic dynamics of test scalar fields. 
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the leading IR effects, and then conhrm that the semiclassical description of the scalar held 
leads to the same equation if we choose an appropriate time coordinate. 

In the latter part of this paper, we remove the constraint of the slow-roll parameter. In 
a general homogeneous, isotropic and spatially hat space, the wave function of a massless 
and minimally coupled scalar held has been derived d- In Section |U by reference to the 
previous study, we conhrm that the propagator has a growing time dependence in a general 
accelerating universe. In Section |5l we extend the nonperturbative analyses in Section |3] in 
a general accelerating universe. We conclude with discussions in Section [6l 


2 Free scalar field theory in an accelerating universe 
with e = 0 

In this paper, we consider a homogeneous, isotropic and spatially hat universe: 

= —dt'^ + a?{t)dy? 

= a^{ri){—dri‘^ + dx^), ( 2 . 1 ) 

where t is a cosmic time and is a conformal time. The dimension of the spacetime is taken 
as D = 4. We consider an expanding era when the Hubble parameter is positive 

H = ->0. (2.2) 

a 

In our notation, the dot denotes the derivative with respect to t. In order to describe 
deformations from dS space, we introduce the slow-roll parameter whose parameter region 
is given by 


e = ^, 0 <£< 1 . ( 2 . 3 ) 

Here the lower limit corresponds to dS space and the upper bound means an accelerating 
universe. 

In Section [2] and [3l we set the slow-roll parameter as it does not change in time 

e = eo : const. (2.4) 

On the background, the scale factor is written by the cosmic time and the conformal time 
as 

+ = ( 2 . 5 ) 
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In the parameter region 02.3 jl . the conformal time approaches zero at the inhnite future. The 
Hubble parameter and the Ricci scalar are given by 

H = Hoa-^\ R = 6{2-e)H^. (2.6) 

Here Ho denotes an initial value of the Hubble parameter. 

Under the scale transformation: 

7] —)■ Ct], X —)■ Ux, (2.7) 

the Hubble parameter and the metric of the background fl2.5p scale as follows 

H C^H, ds^ C~^ds^. (2.8) 

The scaling is understandable since the dimension of each quantity is given as [H] = 1, 
[ds^] = —2. At quantum level, if a vacuum expectation value (vev) of a D-dimensional 
operator [O] = D is written as 

{0{x)) oc (2.9) 

we call that the scaling law is respected. Actually the IR effects which we discuss in Section 
[2] and [3] can be identihed as breakdown of the scaling law. 

We consider a massless and minimally coupled scalar held, the quadratic action of which is 
given by 


^2 = y (2.10) 

For the free scalar held (po, we can expand it by the annihilation and the creation operators 
Op, aj, as 


M^) - j [ap0p(^) + 

(2.11) 

0p(a;) = a-^ x 

(2.12) 

3 — Co 3 Co 

2(l-eo) “ 2 + 1-eo’ 

(2.13) 


where Hul^ is the Hankel function of the hrst kind. The subscript 0 of z/q is inserted to 
emphasize its time independence. We consider the Bunch-Davies vacuum ap|0) = 0 in 
this paper. It should be noted that p denotes the comoving momentum and the physical 
momentum is given by P = p/a. 
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At sub-horizon scales: P ^ H ^ —pr] ^ wave function approaches that in 

Minkowski space up to the conformal factor 


0p(x) 


-1 


X 




• I • • 2 izg -|-1 

^—ip-rj+ip-x—i —^—7 


(2.14) 


In this paper, we focus on contributions from super-horizon scales rather than those from sub¬ 
horizon scales. At the IR scales, specific behaviors for curved spacetimes become dominant. 


In investigating the IR behavior at super-horizon scales: —pi] -C following expansion 

of the Bessel function is useful 

-r{ cos(7ri/o)Vo(j) - (2.15) 

sm(^7rz/oj 


n=0 


(-1)" pY" 
n!r(l ± vq + ny 2 ^ ’ 


(2.16) 


where the double-sign corresponds. The leading IR term of the wave function comes from 
the zeroth order of the expansion of 


X ~ a 


-1 


X 


-v) 


2"''-'r(z/o) 


3+iP'X 


TT 


[-pr]) 




(2.17) 


From fl2.17p . the contribution from super-horizon scales to the propagator is given by 


(<do(^)) - 


22!/o-2p5 




TT 


-v) 


' p<Ha 


d?p 1 
(27r)3 {—pp)‘^‘'o 


(2.18) 


As seen in fl2.13p . uo > 3/2 and therefore the propagator has an IR divergence. The upper 
bound of the momentum integral should be on the horizon since the spectrum is dominant 
at super-horizon scales. To regularize the IR divergence, we introduce an IR cut-off po which 
hxes the minimum value of the comoving momentum. In these settings, the IR contribution 
to the propagator is given by the growing time dependence: 


-2,_U 22‘-»-=>r2(,2„) 


iVoiz:)) zz 


TT-" 


a {-p) 


(2 - 2eo)2"°-3r2(i2o) 


-2 d{-pp) 

pTj^ 21^0—2 

dx 


TT"’ 


'-VOV 


/ N 7 . 21^0 — 2 

'-(1-60)P07? 


(2-2eo)2"°-3r2(i/o) 


TT"’ 


{(l-eo)hfo}V2^° 
1 - eo 


X 


2en 


■{(a/ao)^'° - 1 }. 


(2.19) 
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Here we introduced an initial time Oq which is dehned by the IR cut-off: Oq = {po/Hq)^-^o . 

Physically speaking, a/po means a size of the Universe and it expands with the scale factor. 
Let us recall that super-horizon scales can be represented by the inequality: 

a/pQ > 1/H pq < a. (2.20) 

In accelerating universes a > 0, the degrees of freedom at super-horizon scales increase 
with cosmic evolution. That is one reason why the propagator acquires the growing time 
dependence after the momentum integration fl2.19p . 

We should also note that if the momentum integration is insensitive to the IR cut-off, the 
time dependence of the propagator is given only by the overall factor: 

{(l-eo)L/o}V2^°oci/l (2.21) 

Namely, if so, the propagator respects the scaling law fl2.9p . As seen in fl2.19p . the momentum 
integration is sensitive to the increase of the degrees of freedom at super-horizon scales. The 
scaling law is broken as it is multiplied by the additional time dependence: 

lUS{(a/a„)=‘«-l}. (2.22) 


Here we mention a massless and non-minimally coupled scalar held. Specihcally the addi¬ 
tional term is introduced to the quadratic action 

I (2.23) 

For this scalar held, the index of the Bessel function is deformed as 


3-60 / (3-6o)2 ^6(2 -eo) 

2(1-eo) (1 - ^ 0 )^ ' 


(2.24) 


If ^ is large enough, the scaling law is respected since there is no sensitivity of the IR cut-oh. 
That is why we consider a massless and minimally coupled scalar held in this paper. 


In interacting held theories with massless and minimally coupled scalar helds, physical quan¬ 
tities may acquire growing time dependences through internal propagators. We call them the 
IR ehects in accelerating universes. The IR ehects at each loop level manifest as polynomials 
in a factor growing with time. At late times, the leading IR ehects come from the leading 
powers of the factor at each loop level. This fact indicates that perturbation theories breaks 
down after a large enough cosmic expansion. In the next section, we discuss how to evaluate 
the IR ehects nonperturbatively. 

Before investigating interacting held theories, we show that the propagator fl2.19p can be 
written as the non-local from: 

(^§(x)) c. /Aloga') //«, (2.25) 

^ Jao 
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where H' denotes Such nonlocal representations are also used in P^TT], At the dS 

limit eo —)■ 0, the integration reproduces the logarithmic dependence of the scale factor which 
is found in IM]: 


( 2 - 26 ) 

In fl2.25p . the scaling law is respected for the integrand but it is broken after the integration. 
The nonlocal representation is useful when we count the leading powers of the growing time 
dependence. 


3 Nonperturbative analyses of IR effects in an accel¬ 
erating universe with e = 0 

3.1 Resummation formula of leading IR effects 

In investigating the IR effects nonperturbatively, a rational approach is to resum the leading 
IR effects to all loop orders. The resummation formula in dS space has been derived by 
N. C. Tsamis and R. P. Woodard [8]. In this subsection, we extend the nonperturbative 
investigation in an accelerating universe whose slow-roll parameter is constant. 

In this paper, we consider an interacting field theory whose nonlinear term is given by a 
potential 


S = j [ - ^g>^''d^ipdu(p - V (<p)]. (3.1) 

As a starting point, we consider the Yang-Feldman formalism: 

ip(x) = <fio(x) -i J G^{x,x')-^V{ip{x')), (3.2) 

where is the retarded propagator: 

G^{x,x') = 9{t - t')[ipQ{x),ipo{x')]. (3.3) 

The Yang-Feldman formalism satishes the equation of motion and describes quantum effects 
from all regions of momentum scale. Our purpose is to derive the equation which is valid up 
to the leading IR effects. 
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For the linear term in the right-hand side of fl3.2|) . we extract the leading IR behavior as 


(po{x) ~ LPq{x) 




2"° ^r(z/o)a 

.2"°-^r(no) {(l-6o)i7o}^ 


+ (h.c.)] 

+*P-X + (h-C.)]- 


(3.4) 


Here we introduced the step function which is nonzero at super-horizon scales. The leading 
IR term of the wave function has no time dependence since quantum fluctuations of a massless 
and minimally coupled scalar held are frozen at super-horizon scales. As seen in fl2.17p and 
fl2.19p . the leading IR behavior of the linear term induces the growing time dependences to 
the Wightman functions: 


{ipa(x)ipo(x')), {ipo(x')ipo(x)). 


(3.5) 


In the Schwinger-Keldysh formalism [HE], perturbation theories consist not only of the 
Wightman functions but also of the retarded propagator. To be more precise, each vertex 
integral contains one retarded propagator due to the causality. The nonlinear term in the 
right-hand side of fl3.2p corresponds to the vertex integral. 


Let us extract the leading IR behavior of the retarded propagator in the nonlinear term. In 
evaluating the retarded propagator, we need to know the real and the imaginary parts of the 
wave function. From fl2.15p and fl2.16p . the leading IR terms of the imaginary and the real 
parts come from the zeroth orders of the expansions of and respectively: 




TT 


71 {—pr]Y° 

, .2"o-iF(no){(l-6o)Lfo}"^ , 

= !_ - ^ ^^^-^ 


a ^(—?7) 2 (— 


TT 


P 


^0 


TT 


p, 


2"°+'^oF(z/o)|(i_eo)L7o}^ 


+1 


a-( 3 -^o)jg+ip-x_ 

(3.6) 


The leading IR behavior of the retarded propagator is given by the cross terms between the 
above two parts. The contribution is spatially local in the sense that its spectrum has no 
momentum dependence and then it is proportional to the spatial delta function 


G^{x,x') ~ 


2i/o(l ~ ^o)Ho 


e{t -1') 


d?p 

(27r)3 


[a'-(3-.o) _ ^-(3-.o)]g+ip-(x-x') 


(3 — eo)iLo 


e{t - t') - x'). 


(3.7) 


Substituting the leading IR behavior of the retarded propagator, the nonlinear term is written 
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as 


dt’HQ^a'^°[l-{a'/af-^°]-^V{(p{t’,^)). (3.8) 

The integration of the second term does not contribnte to the leading IR effects since 
{a'/a)^~''° -C 1 at the far past t' t. Neglecting the second term, the vertex integral is 
evalnated as 

j ^/^d'^x' G^{x,x')-^V{lp{x')) j dt' (3.9) 

Note that the time dependence of the integrand can be expressed by H'. In a similar way 
to fl2.25p , it is nsefnl to represent the integral in terms of log a' as 

(3.10) 

As seen in fl3.7p and fl3.10p . the retarded propagator has no IR divergence in itself while the 
vertex integration of it indnces the growing time dependence. 

For example in ip^ theory, two additional propagators appear with the loop level is increased 
by one. One of them is the retarded propagator and the other is the Wightman fnnction. 
From fl2.25p and fl3.10p . the leading IR effects are given by inhnite series of the following 
factor: 

A J d(loga') H'-^ J d(loga") H''\ (3.11) 

where A denotes a conpling constant. The amplitnde of this factor grows with cosmic evoln- 
tion. Therefore, even if the conping constant is small A <C 1, the pertnrbation theory breaks 
down after an enongh time passed: A J d(loga') H'~^ J d{\oga") H”"^ ~ 1. 

Let ns derive the eqnation which describes the leading IR effects nonpertnrbatively. Applying 
the approximations fl3.4p and fl3.10p . the Yang-Feldman formalism is rednced as 

ip{x) = iPq{x) - ^— [ dt'H'-^-^V{(p{t',x.)). (3.12) 

3 - eo y op 

We shonld note that the time dependence of the nonlinear term appears only in the npper 
bonnd of the integral. Differentiating both sides of 03.121) with respect to t, we obtain the 
local equation: 

p{x) = po{x) - ^))- (3-13) 

As discussed below dH, the leading IR term of the wave function is constant and so the 
time dependence of po appears only through the step function which denotes super-horizon 
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scales. As a direct consequence of this fact, the correlation function of (po is proportional to 
the delta function 


((y9o(f,x)(y9o(f',x)) 



{{I-eo)H‘^ay5{Ha-p)5{Ha-H'a') 
2^-°-^r^(no){(l-6o)ifo}^ 

TT 


(2-2eo)2"°-3r2(no) 




(3.14) 


We call this type of fluctuation a white noise. The equation fl3.13^ - 03.14p is known as a 
Langevin equation with a white noise [IS0. 


The Langevin equation can be translated to the equation of the probability density p, which 
is called the Fokker-Planck equation: 


p(tA) = 4>) + 




(3.15) 






Using the probability density, the vev of any field operator is given by 

/ OO 

d(j) p{t,(j))F{(l)), F\ any function. 

■CXD 

Note that (j) is not an operator but a c-number. 


(3.16) 


(3.17) 


In the free field theory U = 0, the solution of the Fokker-Planck equation is given by the 
Gaussian distribution with the growing variance: 


p(f,0) 


1 




exp 




(3.18) 


= A, fdt' = (2 2t„)2-.-2r2(i.„) ^^ _ ^,2 (3 19) 

Jto ^ J ao 

The probability density reproduces the IR behavior of the propagator 02.251) . The IR cut-off 
dependence appears as an initial time. 

In interacting field theories, the most interesting objects are eventual contributions to phys¬ 
ical quantities. In order to evaluate them, we need to solve the Fokker-Planck equation for 

^The standard Langevin equation describes the time evolution of the position of a particle. Here a 
coordinate space is replaced by a field space. 
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a final state. If we know the final state solntion poo, the eventnal contribntion to the vev of 
any held operator can be evalnated by snbstitnting it in the integral fl3.17jl . 


In dS space cq = 0, we can obtain an analytical solntion of the hnal state for an arbitrary 
stable potential. Assnming that an eqnilibrinm state is eventnally established: 


^ Poo( 0 ), 

the Fokker-Planck eqnation fl3.15^ - 113.16p is written as 


0 = 




Svr^ 502 

The hnal state solntion is given by 


Poo(0) + 


1 d 


d 




(3.20) 


(3.21) 


poo{4>) = N ^exp (^ - 1 /( 0 )), 


3if4 


N = 


d(j) exp 


Stt^ 

3H^ 


vm. 


(3.22) 


For example in theory V = the vev of the potential approaches the satnration valne: 


= !§■ 


(3.23) 


The satnration valne is not snppressed by A <C 1 and so it is snrely a nonpertnrbative IR 
ehect. 


On the backgronnd with cq 7 ^ 0, the time independence of the hnal state solntion does not 
hold trne since the Hnbble parameter is time dependent. For a general potential, we cannot 
obtain an analytical solntion of the hnal state. T. Prokopec fonnd that we can obtain it if 
the potential is written as the following form [T 6 ] : 


V{p{x)) = H V{p{x)), p{x) = 


p{x) 


H 


In snch a model, adapting the ansatz that the scaling law is eventnally restored: 

p(f, 0) y Poo(t, 0) = ^Poo(0)) ^ ~ ~H'’ 

the Fokker-Planck eqnation fl3.15p - fl3.16p is written as 

e„A(p„(.A) 4 ) = 

The hnal state solntion is given by 

Poo(0) = iV"^exp ( + ^ 


^0 72 


-vm.. 


(3 — eo)^o 

N = I dcj) exp ( + 4^0^ — 77 ;— \ A V (0)) . 


^0 (3 — eo)^o 


(3.24) 


(3.25) 


(3.26) 


(3.27) 
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Here the tachyonic term appears due to the existence of a nonzero slow-roll parameter. 


As an example, we consider theory which belongs to the solvable case fl3.24p . Since the 
exponent of fl3.27p is not written as a monomial, the evaluation of the integral fl3.17p is 
complicated. Specihcally, the eventual contribution to the vev of the potential is given by 


(■l/(<p(x)))oo = 




^0 


-T 


(3.28) 


The emergence of eo/Az can be explained as follows. Let us recall that the leading IR effects 
are expressed as inhnite series of the factor fl3.1ip . This factor can be rewritten as 


J d(Ahoga') 


d{Xhoga") 


(3.29) 


and the Hubble parameter is a function of Cologa: H = iLg exp (—cq log a). Therefore, if 
eg -C Az, the Hubble parameter can be moved outside of the integral just like in dS space: 

j d{X^\oga') J (i(AMoga") ~ A log^(a/ag). (3.30) 

The higher orders in fl3.28p mean the corrections to this evaluation. 


In an accelerating universe whose slow-roll parameter is constant, we derived the Langevin 
equation from the Yang-Feldman formalism by extracting the leading IR effects. In the 
derivation process, the loop level and the type of Feynman diagram are not specihed. There¬ 
fore, the stochastic equation describes the leading IR effects nonperturbatively. Especially, 
the hnal state solution of the Fokker-Planck equation is useful in evaluating eventual contri¬ 
butions to physical quantities. 


3.2 Semiclassical description of scalar fields 

Before N. C. Tsamis and R. P. Woodard derived the Langevin equation in dS space by 
the resummation of the leading IR effects, A. A. Starobinsky and J. Yokoyama derived the 
same equation from a semiclassical view point [SE]. Here we conhrm that after a slight 
modihcation, the semiclassical approach is consistent with the resummation formula also in 
an accelerating universe whose slow-roll parameter is constant. 

First of all, we make a naive extension of the previous studies The naive investigation 

leads to an inconsistency with the resummation formula. We discuss how to resolve the 
inconsistency in the latter part of this subsection. 
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(3.31) 


We start with the equation of motion: 


d 1 


( 


)V7(X) 


d 


dip 


V{^{x)). 


From the second order differential equation, we extract the equation which describes the 
IR dynamics at super-horizon scales. Since our interest is about physics at late times, we 
neglect the second derivative terms of fl3.3ip naively as 




(3.32) 


Let us recall that in accelerating universes, more degrees of freedom go outside of the horizon 
P = H with cosmic evolution. In other words, super-horizon modes are sourced by sub¬ 
horizon modes. Based on this perspective, we divide the scalar held into super-horizon 
modes ip and sub-horizon modes ip\]Y'. 

ip{x) = ip{x) -b ip\]\{x), (3.33) 


and identify as a source of ip. To be more precise, we set as 


+ (h.,.)]. (3.34) 


(27r 


TT 


p‘ 


and keep the time derivative of it in the equation fl3.32p : 

d d 

^H^^{^{x) + (puv(4} = 


(3.35) 


In the setting fl3.34p . the time derivative of is equal to that of fl3.4p except for the sign: 

^\Jv{x) = -ipo{x). (3.36) 

By moving it from the left-hand side to the right-hand side, we obtain the following Langevin 
equation: 

p{x) = po{x) - ^^V{ip{t,^)). (3.37) 

The noise term is consistent with the result obtained by the resummation formula fl3.13p - 
fl3.14p . In contrast, the coefficient of the nonlinear term is inconsistent with it except in dS 
space eo = 0. As discussed below, the difference of the coefficient originates in the choice of 
the time coordinate the second derivative with respect to which we neglect. 


The consistency with the resummation formula is obtained in the following way. In terms of 
the e-folding number: N = log a, the equation of motion is written as 




1 ^2 




(3.38) 
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Neglecting the second derivatives with respect to N, x rather than with respect to t, x, we 
obtain the equation different from 03.321) : 

(3 - eo)H^-^(p{x) = (3 - eo)H^(p{x) = —^V{(p{x)). (3.39) 

Taking account of the contribution from sub-horizon scales in the same way as 03.34p . 

(3 - eo)H—{^{x) + </9uv(a;)} = —^V(^{x)), (3.40) 

we obtain the Langevin equation which agrees with the result in the previous subsection 

flCT) - 0mi) : 

In an accelerating universe whose slow-roll parameter is constant, we find that the semiclas- 
sical description of the scalar field leads to the same Langevin equation, which is derived 
by the resummation formula. In the derivation process, it is crucial to neglect the second 
derivative with respect to the e-folding number. 

Actually if the slow-roll parameter is time dependent, the e-folding number is no longer an 
appropriate choice to obtain the consistency with the resummation formula. In Subsection 
15.21 we discuss which choice of the time coordinate is consistent with the resummation of 
the leading IR effects in a general accelerating universe. 


4 Free scalar field theory in a general accelerating uni¬ 
verse 


The degrees of freedom at super-horizon scales increase with time as far as the cosmic 
expansion is accelerated. This fact indicates that physical quantities may acquire growing 
time dependences through quantum fluctuations at super-horizon scales even if the slow-roll 
parameter is time dependent e 7 ^ 0. In this section, we investigate the IR effects in a general 
accelerating universe. 


On a general background whose metric is given by fl 2 .ip . the wave function for a massless 
and minimally coupled scalar field has been derived d-. By reference to the previous study, 
the leading IR term of the wave function is evaluated a^l 


(X ~ —i 


.2"‘-ir(z/,) {(l-e,)iL*a^‘}i- 


,+ip-x 


TT 


pv 


(4.1) 


II For simplicity, the reference time U in m is set at the time of horizon crossing ti in this paper. 
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where the subscript * denotes the time of horizon crossing: = p and v is defined in the 

same way as fl2.13p 


2 (1-e) 2 1-e' 


3-e 3 e 

- 7 = 77 + ^ 


(4.2) 


Since each parameter in fl4.ll) is expressed at the time of horizon crossing, the leading IR 
term of the wave function depends on p rather than on t. This fact is consistent with the 
freezing of quantum fluctuations at super-horizon scales. 

From fl4.ip . the IR behavior of the propagator is given by 



(4.3) 


In the second line, the integral variable is changed as p = H'a'. We should note that the 
integral measure induces (1 — e') as follows 


dp = {1 — e')pd{\oga'). 


(4.4) 


Therefore, the propagator is growing with time as far as accelerated expansion continues. 
Furthermore it should be emphasized that the nonlocal quantum effects induced by the 
propagator cannot be subtracted by local counter terms. 

5 Nonperturbative analyses of IR effects in a general 
accelerating universe 

5.1 Resummation formula of leading IR effects 

The growing time dependence of the free propagator fl4.3p indicates that the breakdown 
of perturbation theories takes place in a general accelerating universe. We extend the re¬ 
summation formula in Subsection 13.11 as it is valid even if the slow-roll parameter is time 
dependent. 

To begin with, we extract the leading IR behavior of the free scalar field as 


ipo{x) ~ (po{x) 



2‘^*-W(z/R {(1 - ] i-'* 

i -- - - 



Op -F (h.c.)]. 

(5.1) 
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In a similar way to 03.41) . the time dependence of (po appears only through the step function 
which denotes super-horizon scales. 

As a next step, we extract the leading IR behavior of the retarded propagator. By reference 
to HU, the leading IR terms of the imaginary and real parts of the wave function are given 
by 


0 p(a;) ~ 


— t 




TT 


pv 


TT 


P 


rt* 


ji-.. Jt 


{ / dt' a' ^ + 


3 — e* 


e+*P". (5.2) 


The time dependence of the IR behavior is expressed by the linear combination of 1 and 
dt' a'~^\ independent solutions of the equation of motion at the limit p —?• 0. From the 
above, the retarded propagator is evaluated as 


G^{x,x') ~ [ [ dt" 

J (27r)-^ Ji, 

= —i9{t — t')[ [ dt” a”~^ — x'). 


(5.3) 


In a similar way to fl3.7p . the retarded propagator is proportional to the spatial delta function 
up to the leading IR behavior. 


Substituting fl5.3p . the vertex integral is written as 

-i j G^{x,x')-^V{ip{x')) 


dt' I dt" a” ^ ] ^R((p(F, x)) 


±11 „//-3 


d 


dip 


= — I dt' a'^ [ I dt 


±11 „//—3 


a — 


'•oo a 

,// „//-3 1 O 


dt" a" ]^R(<p(t',x)) 


dt' (a'^ / dt" a" 


// „//—3\ 


/“dt" a"-^- 




(5.4) 


We should note that the amplitude of /^“ dt" a” ^ monotonically decreases with time and 
so /“ dt" a"~^/ /“ dt" a"~^ -C 1 at the far past t' -C t. Therefore, the integration of the 
second term is negligible up to the leading IR effects: 

J ^/—g'd'^x' G^{x,x')-Pv{(p{x')) — ~ J dt' {a'^ 


dt" a" ^)-^V{(p{t',x.)). (5.5) 


Note that the time dependence of the integral appears only in its upper bound. 
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Applying the approximations fl5.1[) and fl5.5p . and then differentiating both sides with respect 
to t, the Yang-Feldman formalism 03.21) is reduced to the Langevin equation: 


^OO 


d 


ip{x) = ifoix) - {a^ I dt' a' ^l/(v 9 (x)), 


where ip^ is given by the white noise due to the freezing of the IR spectrum: 

" r/i _Nrr2 


(99o(i,x)^o(t',x)) = 


(2vr) 


{(1 - e)H‘^a] 5{Ha - p)S{Ha - H'a') 


X 


22"*-2r2(z/,) 


TT 


P 


,2 j /* 


(2-2e)2^-3r2(z/) 






(5.6) 


(5.7) 


The Langevin equation 05.6^ - 05.7p can be translated to the following Fokker-Planck equation: 

1 , r ,-.w d 


p{t,<j)) = ^AH^^p{t,<j)) + {a^ j dt' a' ^)^(^p{t,<j))^V{<j))y 


d(j) 


(5.8) 


A 


(2^^26r2!l!M(i_e)3 

TT^ 


(5.9) 


The coefficient of the second derivative term is positive during the cosmic expansion is 
accelerated. That is why the propagator grows with time. In the free field theory, the 
variance of the Gaussian distribution: dt' A'H'^ reproduces the time dependence of 

the propagator 04.3p . 


On a general background with e 7 ^ 0, we cannot obtain an exact solution of the final state 
even in the model with dimensionless couplings 03.24p . In a similar way to 03.301) . if the time 
variations of 77, e are negligible during the IR effects grow, we can evaluate the final state 
solution just like in dS space: 

Poo(t,0) ~ iV"^exp j d(loga') 

Nr^ dc^ exp ( - y{a^H d{\og a') 

In such a situation, the eventual contribution to the vev of the potential in ip^ theory is 
evaluated as 

A r°° 1 

(y{p>{x)))^r^-{a^H j d{\oga') V. (5.11) 


(5.10) 
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5.2 Semiclassical description of scalar fields 


Here we make the semiclassical derivation of the Langevin equation in a general accelerating 
universe. If we neglect the second derivative with respect to the e-folding number in the same 
way as fl3.39p . the coefficient of its nonlinear term is given by the replacement Cq —)■ e in fl3.4ip . 
However the result of the resummation derivation fl5.6l) shows the different coefficient: 


a 


3 



1 

(3 - ()H 


+ a“ 



(3 - 


( 5 . 12 ) 


From this fact, we find that the e-folding number is not an appropriate choice of time 
coordinate for e 7 ^ 0 . 


To obtain the consistency with the resummation formula, we extend the semiclassical de¬ 
scription of the scalar held in the following way. For a general choice of time coordinate: 


dT = Hdt, 


( 5 . 13 ) 


the temporal part of d’Alembert operator is written as 


A 


d 

3H- = W 
at 


^2 n + 3HH d 
+ 


aT2 


dT 


(5.14) 


In order to compare directly with we choose the time coordinate T as its friction 

coefficient is constant 


H + 3HH 


/io : const. 


The solution of this equation is given by 


n = 




dt' a' \ 


(5.15) 


(5.16) 


Here the boundary condition is set as the amplitude of T monotonically increases with time. 
At the limit e —>■ 0, the time coordinate corresponds to the e-folding number up to the 
constant scale factor: dT = ^^^dN. 

Mo 


Neglecting the second derivatives with respect to T, x from the equation of motion fl3.3ip . 
we obtain the equation which is different from fl3.39p : 

O /‘OO (A O 

= {a^ dt' —(fix) = -—V{(p{x)). (5.17) 

Note that the time independence of po is crucial in the above discussion while the value of 
/io itself does not contribute to the equation fl5.17p . 
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Identifying the contribntion from snb-horizon scales ip\]\ as the sonrce term of the dynamics 
at super-horizon scales (p: 



_i d 

dt'a'~^) —{ip{x) + (Puy{x)] 




(5.18) 


(puv(2:) = 



ap-F(h.c.)], (5.19) 


we obtain the Langevin equation which agrees with the result derived by the resummation 


formula fl5.6h - fl5.7p : 



(5.20) 


Here we made use of the identity between 05.1 p and 05.19^ : (puv = —T’o- 


6 Conclusion 

In accelerating universes, sub-horizon modes with given comoving momenta go outside of 
the horizon with cosmic evolution and then the degrees of freedom at super-horizon scales 
accumulate. As a consequence of the sensitivity to the accumulation, the propagator of a 
massless and minimally coupled scalar field acquires a growing time dependence. In inter¬ 
acting field theories with the scalar field, some physical quantities acquire secular growths 
which spoil perturbative investigations eventually. 

These IR effects are nonlocal in the sense that they are expressed by time integrals of 
functions of TT, e. Therefore, they cannot be subtracted by local counter terms. Especially 
when the slow-roll parameter is time independent, the nonlocal effects can be identified as 
breakdown of the scaling law. 

In order to evaluate the IR effects nonperturbatively, we extended the resummation formula 
of the leading IR effects in a general accelerating universe. Substituting the leading IR 
behaviors of the free held (pQ and the retarded propagator the Yang-Feldman formalism 
is reduced to the Langevin equation with the white noise term fl5.6p - fl5.7p . The difference 
from the Langevin equation in dS space is expressed as the slow-roll parameter dependence 
of the coefficient of each term. 

The white noise term is a direct consequence of the fact that quantum huctuations of a 
massless and minimally coupled scalar held are frozen at super-horizon scales. The nonlinear 
term is expressed as the local operator dV/d(p{p>{x)). That is because up to the leading IR 
ehects, the retarded propagator is spatially local and the time dependence of the vertex 
integral appears only in its upper bound. 
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The Langevin equation can be translated to the Fokker-Planck equation fl5.8jl - fl5.9[) . The hnal 
state solution of the Fokker-Planck equation is necessary to evaluate eventual contributions 
to physical quantities. T. Prokopec found that if the slow-roll parameter is constant and 
the interaction potential contains only dimensionless couplings, we can obtain an analytical 
solution of the hnal state [15]. In a general case, an analytical solution may not always exist. 
In the situation that the time variations of H, e are negligible during the IR effects grow, 
we derived the approximate solution of the hnal state fl5.10[l . 

Furthermore we derived the same Langevin equation by the semiclassical description of the 
scalar held. In the semiclassical approach, we identify the contribution from sub-horizon 
scales as the noise term. As seen in (15.1 p and fl5.19p . we have only to hip the sign of the 
step function variable to obtain the same noise term with the resummation formula. The 
procedure for the noise term is the same as adopted in [6117119HI2]. 

In contrast, we adopt a diherent procedure in eliminating the second derivative terms from 
the equation of motion. If we neglect the second derivatives with respect to t, x in the same 
way as [51EIISHI2], the nonlinear term of the Langevin equation has a diherent coefficient 
with the resummation formula. To resolve the inconsistency, we choose the time coordinate 
T as its friction term is constant in fl5.15p . In terms of T, we can compare directly 
with Neglecting the second derivatives with respect to T, x, the nonlinear term of the 

Langevin equation has the same coefficient with the resummation formula. 

From these facts, we can conclude that as far as we adopt the appropriate time coordinate 
T, the semiclassical description of the scalar held is equivalent to the resummation of the 
leading IR ehects, not only in dS space but also in a general accelerating universe. 
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